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Abstract 

Even though organic conductors have complicated crystalline structure with low symmetry and large unit cell, band 
structure calculations predict multiband quasi-two dimensional electronic structure yielding very simple Fermi 
surface in most cases. Although few puzzling experimental results are observed, data of numerous compounds are 
in agreement with calculations which make them suitable systems for studying magnetic quantum oscillations in 
networks of orbits connected by magnetic breakdown. The state of the art of this problematics is reviewed. To 
cite this article: A. Audouard, J-Y. Fortin, C. R. Physique (2012). 

Resume 

Les conducteurs organiques sous champs magnetiques intenses : systemes modeles pour la physique 
des oscillations quantiques. Bien que les conducteurs organiques presentent des structures cristallines com- 
plexes et de basse symetrie, les calculs predisent une structure electronique multibande quasi-bi dimensionnelle 
conduisant generalement a une surface de Fermi tres simple. Bien que quelques resultats experimentaux deroutants 
soient observes, les donnees de nombreux composes sont en accord avec les calculs, ce qui fait de ces derniers des 
systemes modeles pour la physique des oscillations quantiques dans les reseaux d'orbites couplees par la rupture 
magnetique. L'etat de l'art de cette problematique est passee en revue. Pour citer cet article : A. Audouard, J-Y. 
Fortin, C. R. Physique (2012). 
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1. Introduction 



Quasi-two dimensional (q-2D) organic metals have generally complicated crystalline structures with unit cell 
involving hundreds of atoms and cell parameters as large as few dozens of nm. Nevertheless, band structure 
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Figure 1. (Colour on-line) Fermi surface (FS) of (a) 0-(ET) 4 CoBr4(C 6 H 4 Cl2) [3], (c) /3"-(ET) 4 NH4[Fe(C204)3]-(C3H 7 NO) [4], (d) 
(ET)g [Hg4Cli2 (CeH5Br)2] [5], (e) (TMTSF)2N03, in the temperature range between the anion ordering and the spin density wave 
condensation [6] and (f ) (BEDO)4Ni(CN)4-4CH3CN [7]. These FS achieve networks of hole-type (pink solid linos) and electron-typo 
(green solid lines) orbits. Arrows indicate the quasi-particles path in magnetic field. Labels ct, /3, a, 6, b — a and stand for few 
classical orbits. S and A arc FS pieces corresponding to forbidden orbits, (b) Energy dispersion of 0-(ET)4CoBr4(C6H4Cl2) [3]. 



calculations predict very simple Fermi surface (FS) which allow to view these compounds as model systems 
for quantum oscillations physics. The interplay between atomic arrangement and electronic structure is studied 
in Ref. [1]. Briefly, donor organic molecules such as BEDT-TTF (bis-ethylenedithio-tetrathiafulvalene, further 
abbreviated below as ET), BEDO-TTF (bis-ethylenedioxy-tetrathiafulvalene) or BEDT-TSF (bis-ethylenedithio- 
tetraselenofulvalene) build up conducting planes with various packing types, so called a, /3, k, etc. as reported in 
Ref. [2]. These planes are separated from each other by insulating, generally inorganic, acceptor planes. Each donor 
molecule bears a positive charge controlled by the anion acceptor charge and the stoichiometry. In numerous cases, 
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which are of interest in what follows, the unit cell involves 4 donor molecules with a +1/2 charge yielding two 
holes per unit cell. FS of these compounds, which bears similarity with the 3D alkaline-earth metals, originates 
from one hole-type orbit (labeled Q in the following) with area equal to the first Brillouin zone (FBZ) area, 
that can be approximated as an ellipse in most cases. In the extended zone scheme, Q orbits overlap and gap 
opening is observed at the crossings due to lifting of degeneracy. Depending on e.g. the strength of tt-tv interactions 
between donor molecules, overlap occurs along either one or two directions, yielding networks of orbits liable to 
be connected to each other by magnetic breakdown (MB) in high enough magnetic field. In the first case, we 
deal with a linear chain of coupled 2D hole orbits. Such a network, an example of which is given in Fig. 1(a), 
is an experimental realization of the famous model proposed by Pippard in the early sixties to study MB [8,9]. 
In the second case a 2D network of compensated orbits is observed. This latter network involves two hole-type 
orbits labeled a and b — a in Fig. 1(f) and one electron- type orbit labeled b. Intermediate case is depicted in 
Fig. 1(c) where Q orbits come close together in one direction. In such a case, a large gap opens around the FBZ 
boundary and the network is composed of one hole-type and one electron-type orbit (labeled a) with the same 
cross section. Analogous type of network is also observed in (BEDO-TTF^ReCU-HbO [1,10] and in compounds 
with less trivial FS genesis due to e.g. FBZ folding [5,11] as reported in Fig. 1(d). It can also be observed in the 
Bechgaard salt (TMTSF) 2 N0 3 (where TMTSF stands for tetramethyl-tetraselena-fulvalene). This q-lD metal at 
room temperature becomes q-2D at temperatures below the anion ordering (Tao = 41 K). In-between Tao and 
the spin density wave condensation temperature (Tsdw = 9.4 K), its FS achieves the linear chain of compensated 
orbits displayed in Fig. 1(e) [6]. 

Nevertheless, even not to mention structural phase transitions that can strongly affect the electronic structure, 
experimental data can hardly be reconciled with the calculations of Fig. 1 in few cases. This is mainly due to the 
extreme sensitivity of organic metals to tiny structural details [1] that can be modified by external parameters such 
as temperature or moderate applied pressure. In that respect, high magnetic field-induced quantum oscillations are 
powerful tools for the study of such FS's. Otherwise, in the numerous cases where the calculations reported in Fig. 1 
holds, MB yields oscillatory features, not predicted by the semiclassical models, the quantitative interpretation of 
which is still in progress. 

In Section 2, we report on few examples of puzzling features of the oscillation spectra observed in organic 
metals with predicted FS such as those reported in Figs. 1(c) and (f). The linear chain of coupled orbits, such as 
that displayed in Fig. 1(a), and compensated electron-hole networks, such as reported in Fig. 1(d) and (e), are 
considered in Sections 3 and 4, respectively, from the viewpoint of the quantum oscillations physics in connection 
with MB. 



2. Puzzling oscillations spectra 

Shubnikov-de Haas (SdH) effect, i.e. conductivity oscillations, and de Haas-van Alphen (dHvA) effect, i.e. 
magnetization oscillations, of multiband organic metals are generally studied in the framework of the Lifshits- 
Kosevich model [12,13]. For a single quadratic band r\ of a q-2D system, Landau levels are given by 

E Vtn = Huj v (n + i) +rx cos(feax) (1) 

where n is the Landau level index, uj v = eB cos9/m v is the cyclotron frequency (e = |e| being the electron 
charge), m v is the effective mass, t± is the interlayer energy transfer integral, a± is the distance between conducting 
planes and 9 is the angle between the direction of the magnetic field B and the normal to the conducting layers. 
The oscillatory part of the Kubo conductivity, in the approximation of independent collisions, and magnetization 
are then given by the semiclassical formulae 

_^ - ^ - V" ^ P T - ( 2lxpTl - ^ ( TujJr > j- 27T P k B T D 

P >i 



= 7T~ = Jl ~~E H R ™ cos(2ttp ) (2) 

and 



m osc = -\± " LJ o [-z -R„, P sin(27rp ' ), (3) 

^ Tvprriri \ nw v ) BcosO 

respectively, for SdH and dHvA oscillations, where ctq = e 2 r\a^m ri /(8n 2 h 3 k B T D ) is the zero-field conductivity 
(a detailed calculation is given e.g. in Refs. [14,15,16]), m osc is the oscillatory part of the magnetization [17]. To is 
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the Dingle temperature (To = h/ (2nkBr), where r is the relaxation time). The factor (hLo v /nT± + 2-KpksTD /t±) 
appearing in Eq. 2 is a contribution from the imaginary part of the retarded Green function |SiG 7 , irl (i?)| 2 coming 
from the Kubo formula, where G v ,n{E) = (E — E v , n + iA rltrl (E))~ 1 . The squared term, when expanded, has one 
pole of first order, giving a contribution in l/\A Vin (E)\, and a pole of second order, giving the second contribution 
equal to 2np/(huj v ). In the usual Born approximation the imaginary part |A^, n (_E)| is replaced by its averaged value 
< \A Vtn (E)\ >= nksTo which defines the Dingle temperature. Moreover, the oscillating part of the self-energy, 
which can be included as corrections of 7rfcsTb, is neglected. In the three-dimensional limit where T perp 2> fru) v , 
these oscillations may have significant contribution to the SdH effect. In the clean limit, fkjj v S> ksTo they can 
be discarded. The oscillation frequency F v = fj,m v /(he) of the classical orbit r/, where /x is the chemical potential, 
is proportional to the cross section area A v of the orbit F v = hA v /(2ire), which is physically the quantum flux 
h/e through A v divided by 47T 2 . The damping factor R V , P can be factorized as 

Rv>p{E> B) = Rrj ,pR-q ,pRri .p Rr},p (4) 

where the thermal, Dingle, MB and spin damping factors are given by the expressions: 

,T UoTpiTlri 



V >P 



B cos 9 sinh 



uoTprrin / (B cos 9) 



r,D _ 



uoTr>pm v / (B cos 9) 



(5) 



(6) 



Kp" ^(ipoTHlo)^, (7) 

*- = C0S ^fc2P' ^ 

respectively [9], where Mo = 2ir 2 m e kB /{eh) = 14.694 T/K. In the following, effective masses m n are expressed in 
units of the electron mass m e . Integers and are the number of MB tunnelings and reflections, respectively, 
encountered by the quasiparticle along its closed trajectory, g* is the effective Lande factor. The MB tunneling and 
reflection probabilities are given with a good approximation by po = exp(— Bq/2B) and g 2 = 1 - p 2 , respectively, 
where Bo is the MB field. The exact semiclassical expression for the tunneling parameters can be evaluated by 
considering more generally the Riemann surface of the band structure [18] which is constructed from a polynomial 
in the complex plane C 2 . The singularities of this polynomial are essential in determining the genus of the Riemann 
surface and its homotopy group. Each element in the group is associated to a unique complex quantity S v (which 
is a semiclassical action on this surface) whose real part corresponds to F v and imaginary part to the MB field Bq , 
respectively. Even though direct application of this method is difficult in multiband systems (see few examples in 
Ref. [18]), it gives a nice framework and clear picture of the tunneling process. 

Furthermore, FS warping can lead to beating features in the oscillatory spectra due to the interlayer transfer 
integral tj_ entering the Landau spectrum (see Eq. 1). Warping is accounted for by the warping factors involving 
t± in Eqs. 2 and 3. However, detailed analysis of its contribution to SdH oscillations have demonstrated that t± 
values of a few tenth of a meV, such as it is observed in /3-phase compounds [19], yield oscillatory spectra that 
cannot be explained by the simple addition of a reduction factor [16,20,21,22]. Nevertheless, t± is by one order of 
magnitude smaller in numerous q-2D organic metals, leading to warping factor values close to 1, as far as dHvA 
oscillations are concerned, and/or slowly varying with the magnetic field. 1 

Finally, it should be mentioned that SdH oscillations are generally measured through interlayer magnetoresis- 
tance (R zz ) oscillations. In the case where oscillations amplitude is small compared to the background resistance 
Kb g , it can be assumed that -Aa/ao — AR ZZ /R(, 9 . Besides, contrary to magnetization which, as a thermodynamic 
parameter, is only sensitive to the density of states, magnetoresistance oscillations may originate from quantum 
interference phenomena (QI) [24,25] as early reported for k-(BEDT-TTF) 2 Cu(NCS) 2 [26,27,28]. In such a case, 
Eqs. 2 and 5 to 8 still hold except that the effective mass entering Eqs. 5 and 6 is the difference and sum, re- 
spectively, of the two partial effective masses of each of the interferometer arms. This can lead to small and even 
zero-effective masses for symmetric interferometer. In such a case, oscillations are observed up to high temperature, 
in a range where the oscillations amplitude damping, still accounted for by R^ p , is due to temperature-dependent 
relaxation time governed by inelastic collisions. 



1. For tj_ = 0.04 mcV and = 3.3, which hold for the a orbit of k-(ET) 2 Cu(NCS) 2 [23], Jo(27ttj_ / hu) v ) varies from 0.88 to 
0.98 as the magnetic field increases from 30 to 60 T. 
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Figure 2. (Colour on-line) (a) Oscillatory part of the magnctorcsistancc data at ambient pressure and (b) corresponding Fourier 
spectra of ^"-(ET)4(H30)[Fc(C20*4)3]-C6H4Cl2- As displayed in the inset of (b), SdH oscillations arc still observed at 32 K. Curves 
in main panels of (a) and (b), are shifted down from each other by a constant amount, (c) Temperature dependence of the b 
oscillations amplitude at 30 T and 34.3 T (i.e. 30 T/cos(29°)), respectively, for SdH and dHvA data collected at ambient pressure. 
Solid squares arc SdH data deduced from (a); solid triangles are dHvA data for 6 — 29°. Solid lines are the best fits of the LK 
model assuming either only one closed orbit contributes to the amplitude (dHvA data) or the coexistence of a closed orbit and 
a quantum interferometer with a zero-effective mass and a temperature-dependent scattering rate (SdH data). Each of these two 
contributions arc displayed as dashed lines. From Refs. [33,34]. 



In the following, we will consider a family of organic conductors illustrating the sensitivity of the electronic 
structure, hence of the quantum oscillations spectra, of organic metals to small changes of the atomic structure 
induced by either chemical substitutions or moderate applied pressure. 

Charge transfer salts with the generic formula /3"-(ET)4(A)[M(C204)3]-Solv (where A is a monovalent cation, M 
is a trivalent cation and Solv is a solvent) have raised great interest for many years [29] . These salts which share the 
same /3" packing of the ET molecules can be either orthorhombic, in which case they are insulating, or monoclinic 
q-2D metals. Among these latter salts, denoted as (A, M, Solv) hereafter, many different ground-states, including 
normal metal, charge density wave, superconductivity, and temperature-dependent behaviours can be observed. 
The FS of (NH4, Fe, C3H7NO) is displayed in Fig. 1(c) where the area of the compensated orbits amounts to 
8.8 % of the FBZ area [4]. Analogous FS is also predicted for (H3O, Fe, C6H5CN) [30]. In qualitative agreement 
with these calculations, only one Fourier component with frequency F = 230 T, corresponding to 6 % of the FBZ 
area, is observed in SdH spectra of (H3O, Ga, C6H5NO2) [31]. However, the FS of other compounds of this family 
can be more complicated since four frequencies corresponding to orbits area in the range 1.1 to 8.5 % of the FBZ 
area are observed for (H3O, M, C5H5N) where M = Cr, Ga, Fe [32]. In this latter case, a density wave ground 
state, responsible for the observed strongly non-monotonous temperature dependence of the resistance, has been 
invoked to account for this discrepancy. In contrast, only two frequencies are observed for (H3O, M, C6H5NO2) 
where M = Cr, Ga [35]. Besides, moderate applied pressure have a drastic effect on the SdH oscillations spectra of 
(NH4, Cr, C3H7NO). Namely, whereas up to 6 Fourier components corresponding to orbit area in the range 0.1 % 
to 7 % of the FBZ area are observed at ambient pressure, the spectrum simplifies under moderate pressure since 
only 3 components (Fi = 68 ± 2 T, F2 = 238 ± 4 T, F3 = 313 ± 7 T) corresponding to orbits area in the range 
1.7 % to 7.8 % of the FBZ area are observed at 1 GPa [36]. These scattered results demonstrate the sensitivity of 
the electronic structure of organic metals, hampering any interpretation in the framework of the band structure 
calculations. The last result could nevertheless be interpreted on the basis of the data in Fig. 1. Indeed, as pointed 
out in Ref. [4], the Q orbits of Fig. 1(c) may also intersect along the small axis, leading to one additional orbit 
as reported in Fig. 1(f). In such a case, 3 Fourier components linked by a linear combination settled by the orbits 
compensation, Ff, = F + F(,_ a , should be observed. This picture holds, not only for the above compound under 
pressure for which F3 — Fi + F2 within the error bars, but also for (NH4, Fe, C3H7NO) in the applied pressure 
range from ambient pressure to 1 GPa [37,38] and for (H3O, Fe, C6H4CI2) [33,34,39]. 

Further information can be obtained from the field and temperature dependence of the oscillatory spectra (see 
Eqs. 2 to 8). SdH and dHvA oscillations of the latter compound have been considered as displayed in Fig. 2. 
Whereas a and b — a SdH oscillations follow the LK behaviour, a kink is observed in the mass plot of the b 
oscillations at ~ 7 K. As a result, despite a rather large Dingle temperature (Td = 4 ± 1 K), SdH oscillations 
are still observed at 32 K which, to our knowledge, constitutes a world record for an organic metal. In contrast, 
the dHvA oscillations are in agreement with the LK formula in all the explored temperature range as displayed in 
Fig. 2(c). Keeping in mind that dHvA oscillations are only sensitive to the density of states, it has been shown that 
both the field and temperature dependence of the b oscillations are consistently accounted for by the coexistence 
of a closed b orbit and a quantum interferometer with the same cross section [33,34,39]. Unfortunately, existence 
of a QI path is not consistent with the FS of Fig. 1(f) and this amazing result therefore remains unexplained. 
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F (kT) T (K) 

(a) (b) 

Figure 3. (Color on-line) (a) Fourier analysis in the field range 35-55.3 T of the magnetic torque data, relevant to 
0-(ET)4CoBr4(C 6 H 4 Cl2), reported in the insert. Solid triangles are marks calculated with F a (8 = 0) = 0.944 kT and F p{6 = 0) = 
4.60 kT. Data have been shifted down from each other by a constant amount, (b) Temperature dependence (mass plot) of the 2a 
amplitude. Symbols are data deduced from torque measurements of Fig. 3(a) at 48.3 T. Solid lines are deduced from Eq. 13 with 
m a — 1.81, mp — 3.52, Tjj — 0.79 K and Bq — 35 T. Spin damping factor (see Eq. 8) is either neglected (R^ — 1, p — 1 to 3) or 
accounted for by various g* values. All the data arc normalized to their values at 12 K. 

Otherwise, the frequency F b+a is observed both for (NH 4 , Fe, C3H7NO) [37,38] and (H3O, Fe, C6H4CI2) under 
pressure [34,39]. Taking into account the opposite sign of the a and b orbits, the b+a component cannot correspond 
to a MB orbit [40] in the framework of the FS of Fig. 1(f). In addition, its effective mass is lower than both m a , 
iiii,_ a and m;, which is at odds with the Falicov-Stachowiak model as well. Therefore, this component has not 
a semiclassical origin and can be attributed to the frequency combination phenomenon considered in the next 
sections. 



3. The linear chain of coupled orbits 

MB phenomenon has been intensively studied in 3D alkaline-earth metals over the sixties and seventies (see [9] 
and references therein). In order to compute MB, Pippard introduced the linear chain of coupled orbits in the 
early sixties [8]. The energy spectrum of this model FS, i.e. the density of states in magnetic field in which all the 
allowed classical orbits contribute, can be computed analytically. In contrast to the semiclassical picture, a set of 
Landau bands broadened by coherent MB is obtained instead of discrete Landau levels, due to the traveling of 
quasiparticles on the non-quantized q-lD sheets of the FS. However, the Falicov-Stachowiak semiclassical model 
[40], based on the LK model and involving discrete Landau levels, yielded analytic tool that successfully accounted 
for the field and temperature dependence of the oscillations spectra of these multiband metals [9]. This success 
outshined the Pippard result until the discovery of frequency combinations, 'forbidden' in the framework of the 
Falicov-Stachowiak model, in magnesium [41] and, later on, the experimental realization of the FS proposed by 
Pippard in the starring compound k-(ET)2Cu(NCS)2 [42]. Many other organic metals share this FS topology, an 
example of which is given in Fig. 1(a). The main feature of the dHvA spectra of these compounds is the presence 
of the forbidden Fourier component j3 — a [43,44,45]. Furthermore, depending on the field and temperature 
range explored, effective masses relevant to frequency combinations linked to MB orbits can be at odds with 
the predictions of the Falicov-Stachowiak model, both for dHvA [44,45] and SdH [28,46] oscillations, even taking 
into account QI in the latter case. Numerical and analytical analysis of this spectrum yield, through numerical 
resolution, a non-LK behaviour of the oscillations amplitude [47,48,49]. Besides, the forbidden Fourier components 
observed in these compounds result not only from the formation of Landau bands but are also due to the oscillation 
of the chemical potential enabled by the 2D character of the FS [47,49,50,51,52,53,54]. 



6 



Comprehensive analytical calculations of the oscillations amplitude, taking into account both coherent MB and 
oscillation of the chemical potential, have been reported in Refs. [55,56] and implemented in dHvA data of n- 
(ET)2Cu(NCS)2. However, easy to handle analytic tools necessary to quantitatively account for the experimental 
data were still lacking. Indeed, as pointed out in Ref. [56], the equations are complex and the data could be 
analyzed only numerically. In Ref. [3], a model based on the field-induced chemical potential oscillations, yielding 
analytical non-LK expressions for the Fourier components amplitude is proposed. This model accounts for the 
temperature- and field-dependent quantum oscillations spectra of the strongly two-dimensional charge transfer 
salt ^-(ET^CoB^CetUCb) with an excellent agreement. Remarkably, the unit cell of this compound involves 
two different donor layers [57]. One of them is insulating whereas the FS of the other, displayed in Fig. 1(a), is 
an illustration of the Pippard's model. Fourier spectra of the dHvA oscillations displayed in Fig. 3(a) reveal the 
existence of both classical orbits a, f3, 2/3 — a and forbidden frequencies such as j3 — ot and its harmonics. 

Considering a two band system with effective masses m (i) and band extrema A (i), as displayed in Fig. 1(b), 
the grand potential Q of a 2D slab with area A is given by a series of harmonics p of classical frequencies F v 

^gn ( T, M ) = = Ao) 2 - - Ai) 2 

+ £ £ ^iWT, B) cos(2.p^ + m ) (9) 

p>l v 

in the low temperature limit and within the semiclassical approximation. The chemical potential p and energies 
Aq^j are expressed in Tesla. The universal coefficient $o u o/&b (where $o = h/e is the quantum flux), in unit 
of which the physical quantities are expressed for simplification, is equal to 47r 3 m e /e 2 ~ 4.4 x 10 9 T 2 m 2 /J. r\ 
stands for the fundamental closed orbits that are not harmonics but can be composed of different parts of the 
FS thanks to MB. In particular, two fundamental orbits with frequencies F a = mi(/i - Ai) (i.e. m a — mi) for 
the smallest and Fp = roo(/i - Ao) + mi(/i - Ai), which can be written as Fp = mp(\i - Ap), for the largest 
are defined (see Fig. 1(a)). Within this framework, mo + mi is identified with the effective mass mp of the MB 
orbit p. Each orbit 77 has an effective mass and yields the frequency F n — m v (fi - A v ) where and A^ 
are combinations of the fundamental parameters mo,i and Ao,i. ^17 i s the phase or Maslov index determined 
by the number of turning points on the trajectory r\. At each turning point is associated a value 7r/2, yielding 
^Pa = fp = n f° r orbits 01 and j3 of Fig. 1(a), which is the value for a parabolic band model near the Fermi 
energy. For more complicated orbits, tp v is a multiple of n. For example the phase of the fundamental orbit a + f3 
is ip a +p — 27r. C v is the symmetry factor of the orbit 77, which counts the number of non-equivalent possibilities 
for an orbit to be drawn on the FS. For example, C a = Cp = C2p~ a — 1 and Cp+ a = 2. For a FS composed 
of a single orbit, C v = 1. We notice that frequencies F v depend explicitly on the chemical potential /j. Indeed, a 
frequency physically represents a filling factor proportional to the orbit area in the Brillouin zone or, equivalently, 
the occupation number of the quasiparticles. For example, Fp corresponds to the total area of the FBZ. Eq. 9 
can be easily extended to more complex multiband systems. It is derived from the usual semi-classical technique 
using the Poisson formula applied in the case where hui c is small compared to the chemical potential fi [17]. The 
oscillating term proportional to B 2 in Eq. 9 contains all the possible contributions of MB between the bands. In the 
Grand Canonical Ensemble, the chemical potential is independent of the magnetic field, leading to the expression 
of the magnetization m = — dCl/d(B cos 9), hence Eq. 3. In the case where the electron density is fixed, which 
is common in dHvA experiments, the chemical potential generally depends on the magnetic field and oscillates. 
These oscillations can be damped if the system is connected to a reservoir of uniform electron density [47]. The 
electron density per surface area n e is defined by dfi/d/j = -n e . In zero-field, Eq. 9 yields n e = (mo + mi)/jo - 
moAo - miAi where /jo is the zero-field chemical potential. In the presence of a magnetic field, /i satisfies instead 
the following implicit equation 

F = W-EE ^^C V R V , P (T, B) sin(-g^L + m) . (10) 
— mp-Kp Bcoso 



The small oscillating part of the magnetization can be computed systematically by inserting Eq. 10 in the 
/n-dependent terms of Eq. 9, in particular the frequencies, and computing m = — dCl/dB cos (9). A numerical reso- 
lution of this resulting implicit equation could be done to obtain recursively the field dependence of fi. Nevertheless, 
a more user-friendly controlled expansion in powers of the damping factors R V , P can be derived systematically at 
any possible order. Up to the second order, the following analytical expression is obtained for the oscillating part 
of the magnetization: 
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EE^,(WS+^)+E E ^w^cox 



Trprriri \ B cos # / z — ' z — ' ivp'mg 

r] p>l T],7]'p,p'>l 

P F V + P 'F V , , \ . ( pF v -p'F v , 



(11) 



where frequencies F v = m v (fJ,o — A 7) ) are evaluated at /i = ^o- Eq. 11 can be written as a sum of periodic 
functions m oac = 'Y^ i A i sm(2-nF i /(B cos9)) where the index i stands for either classical orbits r\ or forbidden 
orbits such as j3 - a. For the most relevant amplitudes, from the experimental data viewpoint, we obtain the 
following expressions: 

F a F a F a 2F a 

A a = Ra,l + Ra,lRa,2 + Ra.2Ra.:i H Rb,xRix+B,1, (12) 

nm a 2nmg onmg nmg 

A20, — — — Ra,2 H —Ra 1 — 7; —Ra,lRa,3, (13) 

2nm a irmp dirmg 

As = -^-Rb,i + t^—R0,iRb,2 + -^—Rb,2Rb,s + ^-Rb,iRc+b,i, (14) 

A28 = ~^-R B ,2 + —R 2 B,1 ~ ^-^R B ,lR B ,3, (15) 
Zirmg Jimp Sttttib 

Ag-u = — —Ra.iRgA ^—^-Ra.2R a +B,i — —Rg,2Ra+g ,1 , (16) 

irrng Tvmg imig 

Ag+a = — - — Rg+a.i H — ^—Ra,iRg,i, (17) 

■Kmg +a Tvmg 

A2g~u = — 2;3 - R2B-a,i + 2/3 " R a ,iRg,2 + -rr^——Ra,sRa+B,2- (18) 
irm2B-a ZirmB vnrriB 

These equations differ from the LK model in the sense that a basic orbit such as a, which should involve only one 
damping factor A a = F a R a /{Trm a ) according to Eq. 3, involves additional terms which are power combination of 
different damping factors. However, deviations from the LK behaviour due to the high order terms are significant 
only in the low temperature and high field ranges. This statement also stands for the classical orbit 2/3- a since 
R2g- a is significantly higher than the product R a ,iRg,2 (see Eq. 18). Eq. 11 bring out forbidden frequencies such 
as P — a (see Eq. 16). The amplitude of such Fourier component arises from the combinations of classical orbits 
a and /3, hence only at the second order, through the damping factor product R a ,iRg,i- More generally, the 
main rule that can be derived is that products of damping factors involved at a given order represent algebraic 
combinations of frequencies. For example, the orbit /3 in Eq. 14 can be viewed as the combinations 2/3 - /3, or 3/3 - 
2/3, yielding the factors Rg^Rg,i and Rg^Rp^, respectively, as additional factors entering the amplitude Ag. As a 
consequence, forbidden frequencies arise from the order two since they cannot be due to any single classical orbit. 
Algebraic sums of damping terms in Eqs. 13 and 17 may possess minus signs which account for n dephasings, and 
may cancel at field and temperature values, depending strongly on the effective masses, Dingle temperature, MB 
field, etc. as displayed in Fig. 3(b) relevant to 2a. This point has already been reported for the second harmonic 
of the basic orbits both for compensated [58] and un-compensated [54] metals. 

Finally, even taking into account the contribution of QI, stronger deviations from the semiclassical model are 
reported in the case of SdH spectra of several compounds with this FS topology [3,47,46]. These features are not 
accounted for by the above calculations which therefore requires a specific model. 



4. Networks of compensated orbits 



Let us consider now the other class of FS's, widely encountered in organic metals, which are built up with 
compensated electron and hole orbits (few examples are reported in Figs. l(c)-(f)). Depending on MB gaps 
value, either isolated orbits, ID or 2D networks are observed. As mentioned in Section 2, few hints of forbidden 
frequency combinations are observed in SdH spectra of /3"-(ET)4NH4[Fe(C2C>4)3]-C3H7NO. Besides, QI which 
requires MB gaps crossing as well, have been tentatively inferred to account for the oscillatory data of /3"- 
(ET)4(H30)[Fe(C204)3]-C6H4Ci2 reported in Fig 2. However, to our knowledge, the only compounds studied 
from this viewpoint belong to the family (ET^pE^X^CeHsY^] (where X, Y = CI, Br) [5,59,60,61,62,63]. 
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Figure 4. (Colour on-line) Forbidden (left side) and MB-induced (right side) S orbit relevant to the Fermi surface of 
(ET)g [Hg4Cli2 (CgHsBr^] [5]. Numbers on the right side indicates the successive steps of the quasiparticles starting from an 
(arbitrary) origin marked by the open circle. Arrows mark the quasiparticles path. Pink and green arrows correspond to hole and 
electron parts, respectively, of the Fermi surface. Black arrows and circles mark the magnetic breakdown tunneling junctions and 
reflections, respectively. Taking into account the opposite sign of the compensated electron and hole basic orbits (labeled a), the 
considered MB orbit yields a frequency corresponding to the area of 5. 

Indeed, due to moderate MB gaps, compounds with X = CI achieves 2D networks as reported in Fig. 1(d). In 
addition, small scattering rate (To = 0.2 ± 0.2 K) is observed [61]. In line with the Falicov-Stachowiak model, 
Fourier components are linear combinations of the frequencies linked to the compensated closed orbits a and to 
the 'forbidden' orbits 8 and A [5,60,61,62]. This is the case of e.g. 2a + S which is observed both in SdH and 
dHvA spectra. Regarding magnetoresistance, components linked to QI such as a + S and 2a + 6 + A are also 
observed in the spectra. Strictly speaking, there is no forbidden orbit in a 2D network of compensated orbits. 
For example, keeping in mind that electron and hole orbits have opposite signs, a MB-induced orbit with the 
frequency corresponding to 5 can be defined for the FS of Fig. 1(d) as displayed in Fig. 4. However, according 
to the Falicov-Stachowiak model, its amplitude should be very small. Indeed, this orbit involves 6 MB tunnelings 
and 10 reflections leading to small R MB values (see Eq. 7). Besides, its effective mass mg should amount to im a 
leading to small R D and R T values, as well. In violation of these statements, 'forbidden' orbits such as 8 or 4a + 5, 
with a small effective mass (mg ~ 0.4m a , m,4 a+ s ~ 0.3m a [61]) and a large amplitude, are observed. As a result, 
these two Fourier components which should not be observed in the framework of the semiclassical model, are 
preponderant above few kelvins at moderate fields (see Fig. 5). 

Contrary to magnetoresistance data, the whole dHvA spectra is quantitatively accounted for by the Falicov- 
Stachowiak model in the field range up to 28 T [62]. This latter point is at variance with the data for the linear chain 
of coupled orbits considered in Section 3. In order to solve this discrepancy, a linear chain of compensated orbits, 
i.e. the ID network which accounts for the FS in Fig. 1(e), has been considered as a first step in Refs. [58,64]. 
Indeed, its energy spectrum can be easily deduced from the Pippard's method [8]. In this model the effective 
masses linked to the electron m e and hole band are taken to be independent. For simplicity, two quadratic 
potentials have been considered to model the FS, i.e. quadratic functions of the quasi-momenta are assumed for 
the bands dispersion, as within the free-electron model in two dimensions. The bottom of the electron band is 
set at zero energy while the top of the hole band (inverse quadratic potential) is at A > 0, with the possibility 
for the quasiparticle to tunnel through the gap between two successive orbits by MB, i. e. between the pink and 
green parts of the FS in Fig. 1(e). Given an energy E, the k-space areas of the closed electron and hole orbits 
are respectively given by S e = 2-Km R E and Sh = 27rm^(A — E), which are both quantized. The zero field Fermi 
energy Ef is given by the condition of compensation S b (Ef) = Sh(Ef) or Ef = [mh/(m e + mh)]A. The unique 
fundamental frequency of this system is therefore equal to F a — S £ (Ef) /2n = Sh(Ef) /2tv = m e m,hA/(m e + m/i). 

Strikingly, the field-induced chemical potential oscillation, calculated by extremizing the free energy is strongly 
damped compared to the uncompensated case of Section 3 [58,64]. For this reason, Fourier amplitudes can be 
calculated within the LK formalism. Nevertheless, it is important to notice that, in systems with compensated 
bands, there is an infinite number of classical orbits contributing to any given Fourier component since a semi- 
classical trajectory around a successive electron and hole pockets has a zero area. Here orbits with zero area 
can be classified by their increasing masses (I — l)m e + (I — l)mh, where Z > 1 is an integer. The amplitude 
A a of the Fourier component with frequency F a is therefore not only dependent on the orbits composed of one 
electron or one hole orbit with effective mass m e or mu, respectively, but also on a series of additional multiple 
orbits composed of Z — 1 electron and Z — 1 hole orbits plus one electron or one hole orbit, with effective mass 
m e (l) = lm e + (I — l)m,h or m,h(l) = (I — l)m e + Irrih, respectively. A a can be computed exactly in this model, 
within the LK theory, and is given by the expression 
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Figure 5. (Colour on-line) (a) Oscillatory part of the magnetoresistancc and (b) corresponding Fourier spectra of 
(ET)g [Hg4Cli2 (C6H5 Br)2] at various temperatures [61]. Fourier spectra of (c) magnetoresistance and (d) do Haas-van Alphcn 
oscillations for various directions of the magnetic field [62] . Data in (c) and (d) have been measured simultaneously. Labels a and 
S correspond to closed and forbidden orbits, respectively, depicted in Fig 1(d). Contrary to magnetoresistancc oscillations spectra 
which exhibits forbidden frequency combinations, all the observed dHvA components arc accounted for by the Falicov-Stachowiak 
scmiclassical model. 
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Figure 6. (Color on-line) (a) Fourier spectra of the dHvA oscillations, displayed in the insert, of a linear chain of compensated orbits, 
for various values of the magnetic breakdown field Bq. The data arc calculated from numerical resolution of the field-dependent 
free energy equation with A — 1 and effective masses m e — 1 and — 2.5 (see text and Rcf. [58]). F a — 7/2 is the fundamental 
frequency for this model, (b) B/T dependence of the first harmonic amplitude A a for various values of the tunneling probability 
(hence of the magnetic field). Solid symbols come from numerical resolution of the scmiclassical spectrum, which fall onto the solid 
lines corresponding to the Lifshits-Koscvich approximation given by Eq. 19 including contributions of zero area orbits. Dashed lines 
correspond to the first order term / — 1 in Eq. 19 or cquivalently to Eq. 3. 
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where the combinatorial quantities can be defined by the summations 
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(20) 



These positive integers count the number of non-equivalent orbits for a given mass m e ^(l) visiting 2n pockets 
on the path. Table 1 displays the numerical values of S(l, n) computed directly from Eq. 20. Damping factors are 
given by ^(m^jj) = RfjR^Rf where i = e,h. po and qo are the MB tunneling and reflection probabilities (see 
Eqs. 5 to 8). Even though the contribution of large masses is usually negligible at low field or high temperature, 
it can be significant in the large B/T limit in which case they have to be taken into account, provided the Dingle 
temperature is small [58]. Indeed, the temperature dependent factors are close to unity in this limit. 



Table 1 

First values of the coefficients <S(Z, n) representing the number of non-equivalent orbits for a given mass m e ^j (I) with 2n magnetic 
breakdowns, 1 < n < 2(1 — 1). 
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Technically, the summation over all the orbits with zero area is evaluated using an analogy with an Ising model 
on a linear chain in magnetic field. Indeed, it is useful to view each quasiparticle traveling on the linear chain 
represented by the FS of Fig. 1(e) as a particle performing a one dimensional random walk. A set of positions on 
the chain {xi}i = o,2n can be defined for a given periodic orbit, with boundary conditions xo — X2n = 0. These 
coordinates take integer values (negative and positive) and define the pocket inside which the quasiparticle is 
located. In what follows, Xi with i even (odd) is the position of a quasiparticle in one electron (hole) band. A 
closed path has an even number of steps 2n. Given Xi, the particle can also orbit a number rii > of times around 
the surface before going to the next band by MB. We can rewrite coordinates Xi with mean of forward/backward 
variables en = ±1 such as Xi = Xi-i + <Ti(xi-i — Xi-2). Here a = 1 when the particle is going forward and 
di — —1 when it is going backward. An adequate set of variables is given by yi — Xi — Xi-i = ±1 which satisfy 
Vi = (TiUi-i or <7j = y%y%-i. All the possible orbits contributing with a given effective mass can then be counted 
by summing over all the possible yi — ±1 restricted to the boundary conditions. It can then be proved that this 
combinatorial number is equal to the partition function of an Ising model in a (complex) magnetic field and which 
is solvable, leading to Eq. (20). In two-dimensions, compensated networks such as the FS reported in Fig. 1(d) 
could also lead to non-negligible contribution of zero-area orbits, but their explicit evaluation is difficult since a 
similar analogy with the previous calculation would lead equivalently to the computation of the partition function 
for a two-dimensional Ising model in magnetic field. 

The relevance of Eq. 19 is evidenced in Fig. 6(a). Oscillations data in this figure are calculated from numerical 
resolution of the field-dependent free energy equation given in Ref. [58]. As expected, Fourier analysis exhibits only 
one frequency F a and harmonics. Fig. 6(b) compares the temperature dependence of A a deduced from numerical 
resolution (solid symbols) to the predictions of Eq. 3, which neglects any contribution of orbits with / > 1 (dashed 
lines). A growing discrepancy is observed as the temperature decreases, except in the absence of MB (qo = 1) in 
which case only the basic closed orbits contributes. In contrast, an excellent agreement with Eq. 19 (solid lines) 
is observed. Shortly speaking, at variance with uncompensated orbits networks, the LK formula accounts for the 
quantum oscillations in a linear chain of compensated orbits (one-dimensional network), provided the contribution 
of all MB orbits is taken into account for each frequency. 
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5. Conclusion 



Although few puzzling results are observed, as reported in Section 2, numerous organic conductors achieve 
Fermi surfaces that can be regarded as model systems for the study of quantum oscillations in networks of orbits 
coupled by magnetic breakdown. Two types of networks can be distinguished, namely the linear chain of coupled 
orbits and the networks of compensated orbits (see few examples in Fig. 1). These networks are considered in 
Sections 3 and 4, respectively. 

A user friendly analytical tool, taking into account the chemical potential oscillations in magnetic field, has 
recently been proposed to account for dHvA oscillation spectra of the linear chain of coupled orbits [3]. The 
main feature of the derived formulae is, besides standard Lifshits-Kosevich terms, the presence of second order 
corrections liable to account for the observed Fourier components that are 'forbidden' in the framework of the 
semiclassical Falicov-Stachowiak model. This model has been successfully implemented in the dHvA oscillations 
data of ^-(ET^CoBr^CeHUC^). Obviously, this encouraging result needs to be confirmed in the case of other 
relevant compounds. In addition, even stronger deviations are observed for SdH spectra which are still unexplained. 

Regarding compensated orbits, either one- (see Fig. 1(e)) or two-dimensional (see Fig. 1(d), (f)) networks can 
be observed. At variance with the above case, oscillations of the chemical potential are strongly damped in one- 
dimensional networks and the Lifshits-Kosevich theory applies in this case, provided multiple orbits are taken 
into account. Unfortunately, to our knowledge, no experimental data relevant to one-dimensional networks are 
available to check this result. Nevertheless, dHvA oscillations spectra of the two-dimensional network achieved by 
the (ET)g [Hg4Cli2 (CeHsBr^] compound are nicely accounted for by the Falicov-Stachowiak model in magnetic 
fields below 28 T which is in line with the above reported damping of the chemical potential oscillations. However, 
SdH oscillations evidences strong deviations from the Lifshits-Kosevich model, including the presence of 'forbidden' 
Fourier components with high amplitude, besides the contribution of quantum interference. These results still 
require interpretation. 

Regarding new problematics, promising perspective are brought about by organic metals with two different 
donors planes, few examples of which can be found in the literature. In few cases, a metal-insulator transition is 
observed as the temperature decreases [65,66] suggesting two insulating planes at low temperature. Compounds 
with one metallic and one insulating layer with either compensated orbits [67] or, as reported in Section 3, linear 
chain of coupled orbits [3,57] have also been reported. In these cases, a strongly two-dimensional behaviour is 
observed. More appealing case is provided by unit cell with two different metallic layers [68,69,70], hence with two 
different Fermi surfaces [69]. Such a structure could be relevant to the bilayer splitting phenomenon reported for 
cuprate superconductors [71] and addresses the question of a three-dimensional Fermi surface definition. To our 
knowledge, quantum oscillations data at high field have only been reported for #-(BEDT-TSF)4HgBr4(C6H5Cl) 
for which complicated spectra, with no clear frequency combinations are observed [69]. In any case, the oscillatory 
data cannot be interpreted on the basis of a mere addition of the contributions of the two Fermi surfaces pertaining 
to each of the two donors planes. Synthesis of other compounds and further experiments are needed to get a better 
insight on this issue. 
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